// Programa de scilab escrito 23/11/2010
// Como trabalho 8 da matéria FFI 0321 - Modelagem Matemático-Computacional
// no IFSC-USP
// _o_o_ oOo _o_o_


// A) Queda livre, fixando quando encontra adjacente
TAM=10;
matriz = zeros(TAM*2,TAM);


foo=1;
while foo == 1,
  x=int(rand()*TAM)+1 // escolhendo de onde soltar a partícula
  if x>1 then
    esq = -1;
  else
    esq = 0;
  end
  
  if x<TAM then
    dire = 1;
  else
    dire = 0;
  end
  
  esq, dire;
  
  y=1;
  barr=1;
  while barr == 1, // loop para a particula cair
    if (y==TAM*2) | (sum(matriz(y+1,x+esq:x+dire))>0) then
      barr=0;
      y;
    else
      y=y+1;
    end
  end
  matriz(y,x)=1;
  
  if y == 1 then
    foo = 0;
  end
end

// FAZER O PLOT DA matriz KRAIO TTM

foo=meshgrid(TAM*2:-1:1,1:TAM)'.*matriz

//plot(1:TAM,foo')